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[2]  ABSTRACT:  Optical  storage  systems  offer  the  potential  for  drastically  increased  data  transfer  rates 
through  the  use  of  parallel  access.  The  fidelity  and  capacity  of  these  optical  systems  are  impacted  by  the 
interfaces  used  to  access  the  parallel  data.  In  particular,  the  coding  and  signal  processing  components 
of  the  system  must  be  designed  with  the  2D  data  format  in  mind.  This  research  has  focused  on  novel 
methods  of  coding  both  the  object  arms  (e.g.,  interleaving)  and  the  reference  arms  (e.g.,  apodizing)  of 
holographic  optical  storage  systems.  Our  research  efforts  have  produced  several  significant  results:  (1)  An 
information-theoretic  methodology  has  been  established  for  quantifying  the  performance  of  VH  systems, 

(2)  A  multidimensional  interleaving  strategy  has  been  developed  for  optimizing  bit  positions  within  a  2D 
channel,  (3)  An  information-based  design  paradigm  has  been  implemented  for  the  design  of  data-bearing 
optical  systems,  and  (4)  Tradeoffs  concerning  the  balance  of  noise  reduction  with  storage  density  have  been 
produced  for  both  apodized  and  unapodized  systems.  Item  (4)  concerns  work  that  is  not  yet  published  and 
a  copy  of  an  internal  report  describing  this  work  is  attached. 

[3]  TECHNICAL  PROJECT  SUMMARY:  An  outline  of  the  research  tasks  completed  during  the  period  of 
this  contract  is  given  below.  The  detailed  results  associated  with  each  of  these  tasks  have  been  reported  in 
numerous  journal  publications  and  student  theses.  A  comprehensive  list  of  these  is  also  provided  below. 

1.  Information  theoretic  capacity /density 

(a)  identified  noise  and  error  sources  relevant  to  VHM. 

(b)  Discovered  balance  between  number  of  pages  and  information  per  page. 

(c)  Extended  scalar  information  formalism  to  vector  approximation. 

(d)  Computed  VHM  capacity  bounds  for  various  random  and  systematic  errors. 

2.  2D  interleaving 

(a)  Formulated  optimal  interleaver  in  terms  of  topological  framework. 


(b)  Combined  various  codes  with  optimal  and  near-optimal  interleaving. 

(c)  Evaluated  BER  and  capacity  of  interleaved  systems  over  wide  range  of  system  parameters. 

(d)  Experimental  demonstration  of  interleaving  and  ECC. 

3.  Information-based  design 

(a)  Extended  information-theoretic  formalism  to  optical  propagation  and  design. 

(b)  Comptued  the  information  capacity  of  Fresnel  propagation  for  various  detection  schemes. 

(c)  Computed  the  information  capacity  of  coherent  imaging. 

(d)  Investigated  the  optimization  of  pagewise  mutual  information  and  defined  the  information-theoretic 
space-bandwidth-product. 

(e)  Defined  new  optimization  metric  for  information-sensitive  optical  design. 

(f)  Demonstrated  utility  of  new  design  metric  on  various  systems. 

4.  Apodization  tradeoffs 

(a)  Derived  angular  selectivity  profile  in  terms  of  fundamental  system  parameters. 

(b)  Studied  relationship  between  apodization  profile  and  crosstalk  noise. 

(c)  Derived  the  M/#  of  an  apodized  VHM  system. 

(d)  Experimentally  demonstrated  use  of  various  apodization  profiles. 
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CHAPTER  3 
Apodization 


3.1  Introduction 

The  key  to  achieving  an  enormous  storage  capacity  in  a  volume  holographic 
memory  is  the  ability  to  multiplex  a  large  number  of  holograms  within  the  same 
physical  volume.  The  choice  of  a  multiplexing  scheme  dictates  much  of  the  system 
architecture  necessary  to  allow  extraction  of  the  desired  hologram  from  the  col¬ 
lection.  There  are  a  large  number  of  proposed  multiplexing  architectures,19’37-41 
each  of  which  has  its  own  unique  advantages  and  disadvantages.  We  will  pri¬ 
marily  consider  angular  multiplexing  but  much  of  what  is  presented  here  can  be 
easily  applied  to  a  variety  of  other  multiplexing  approaches. 

A  hologram  is  created  by  interfering  the  data-bearing  object  beam  with  a 
reference  beam  at  a  specific  choice  of  the  multiplexing  parameter,  say  ip0.  Later, 
when  we  wish  to  reconstruct  the  desired  object  beam,  we  use  the  reference  beam 
corresponding  to  tp0  resulting  in  a  high-fidelity  reconstruction  of  the  object  beam. 
As  we  perturb  the  multiplexing  parameter,  in  general,  the  irradiance  and  quality 
of  the  reconstructed  image  decay.  We  wish  to  investigate  in  considerable  detail 
precisely  how  the  diffraction  efficiency  changes  as  we  adjust  the  multiplexing 
parameter  away  from  tpo- 

When  there  are  multiple  holograms  stored  in  the  same  memory,  illumination 
with  the  desired  reference  beam  reconstructs  the  desired  page  but  it  also  pro¬ 
duces  contributions  from  all  the  other  holograms  albeit  at  very  low  diffraction 
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efficiencies  relative  to  the  desired  page.  Although  with  proper  selection  of  the 
multiplexing  spacing  the  hologram-to-hologram  cross  talk  can  be  significantly  re¬ 
duced,  it  remains  as  a  noise  source  that  scales  with  the  number  of  participating 
holograms.  In  the  case  of  very  high  density  storage,  we  expect  that  the  cross  talk 
level  will  limit  acceptable  retrieval  of  the  desired  data  page.  It  is  in  this  cross 
talk  limited  case  that  we  wish  to  consider  in  this  chapter. 

We  begin  by  reviewing  Bragg  diffraction  for  thick  media  and  proceed  to  de¬ 
velop  expressions  for  the  diffraction  efficiency  and  hologram-to-hologram  cross 
talk  in  terms  of  the  underlying  system  parameters  specifically  for  angular  multi¬ 
plexing.  The  noise-to-signal  (NSR)  metric  is  introduced  in  sec.  3.5  to  character¬ 
ize  the  level  of  cross  talk  present.  Reference  beam  apodization  during  recording 
and/or  reconstruction  is  described  as  a  technique  to  improve  the  NSR  and  im¬ 
prove  both  the  multiplexing  density  of  the  memory  and  the  total  number  of 
multiplexed  holograms.  In  the  last  part  of  this  chapter  the  relationship  between 
apodization  and  M/#  is  derived. 


3.2  Bragg  diffraction 

The  first  step  in  understanding  volume  diffraction  of  the  reference  beam  is  to 
consider  how  a  single  location  in  the  material  responds  to  the  incident  optical 
wave.  The  first-order  Born  approximation7,42  allows  us  to  characterize  the  re¬ 
sponse  at  each  point  as  a  spherical  scattering  center  with  a  scattering  coefficient 
proportional  to  the  local  change  in  the  material  permittivity  and  amplitude  of 
the  illuminating  wave.  The  field  at  a  point  in  the  output  plane  can  then  be 
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represented  as  the  coherent  superposition  of  the  scattered  wavelets  from  all  the 
contributing  points  in  the  volume. 

Note  that  the  Born  approximation  is  valid  only  for  weak  holograms.  If  the 
diffraction  efficiency  is  significantly  large,  then  the  more  rigorous  coupled  wave 
theory8  must  be  applied.  The  small  diffraction  efficiency  assumption  can  be 
thought  of  in  two  ways.  First,  that  scattering  from  one  location  does  not  scatter 
again  at  some  other  location  and  two,  the  reference  beam  is  not  significantly 
attenuated  due  to  the  coupling  of  energy  into  the  diffracted  beam,  i.e.,  the  energy 
lost  in  the  reference  beam  due  to  the  scattering  is  negligible. 

By  considering  the  reciprocal  wavevector  space  (angular  spectrum),  we  can 
represent  the  superposition  of  scattering  sources  in  terms  of  a  transfer  function 
A(ki,  kd)  that  provides  the  amplitude  of  the  diffracted  field  with  wavevector  kd  due 
to  the  component  of  the  reference  beam  with  incident  wavevector  fcj.  A(kl,  kd)  es¬ 
tablishes  a  convenient  mapping  between  input  wavevectors  and  output  diffracted 
wavevectors.  The  derivation  can  be  found  in  refs.  [42  and  43].  In  the  derivation 
it  is  assumed  that  the  spatially  varying  material  permittivity  can  be  expressed 
as: 


e(r)  =  e0  +  l{r)eiR*f  (3.1) 

where  two  conditions  must  be  satisfied.  First,  the  amplitude  of  the  modu¬ 
lated  permittivity  must  be  much  smaller  then  the  bulk  average  permittivity: 
|e(f)|  <C  eo-  The  second  requirement  is  that  the  spatial  bandwidth  of  the  permit¬ 
tivity  modulation  be  much  less  than  the  grating  frequency.  Both  conditions  will 
be  in  practice  satisfied  for  photorefractive  holographic  storage. 
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Accepting  the  previous  assumptions,  the  wavevector  transfer  function  can  be 
written  as:42,43 


At?  T  \  f 

kd)  —  I  - 
Jv  2iJ> 


(r)  'fEj 


L-2  _  JU2  _  1.2 

where  ^(r)  is  the  amplitude  of  the  reference  field  at  the 
scattering  volume,  V. 


dr  (3.2) 

position  r  inside  the 


Note  that  if  the  reference  beam  only  consists  of  one  wavevector,  fcj,  that 
A(ki,  kd)  defines  the  angular  spectrum  (in  terms  of  the  variable  kd)  of  the  diffracted 
wave.  We  can  then  write  the  diffracted  field  as  the  inverse  2D  Fourier  transform 
of  the  wave’s  angular  spectrum: 

Eout(r)  =  J  J  A(ki,kd)elkdfdkd,xdkd>y  (3.3) 


3.2.1  Ideal  volume  phase  grating 


Understanding  diffraction  from  the  ideal  thick  phase  grating  will  serve  as  a 
starting  point  for  the  later  more  complicated  developments.  Consider  the  simple 
case  where  the  permittivity  modulation  consists  of  an  undamped  single  grating: 
e(r)  =  1  and  the  grating  vector,  Kg,  is  defined  through  the  interference  of  an 
object  monochromatic  wave  with  wavevector  k0  and  a  recording  reference  wave 
with  wavevector  kT: 


Kg  —  k0  kr  (3.4) 

-  If  a  crystal  of  dimensions  WxHxL  (refer  to  fig.  3.1)  with  negligible  absorption 
(an  ideal  assumption)  is  illuminated  with  a  reconstruction  reference  beam  of 
wavevector  we  wish  to  compute  the  amplitude  of  the  diffracted  field  with 
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wavevector  kd.  By  eq.  3.2,  we  know  the  amplitude  of  the  diffracted  field  in  the 
direction  defined  by  kd  is  A{kukd)  which  can  be  simplified  using  the  assumed 


£  =  (kd-k0)  +  (kr-ki)  (3.10) 

To  satisfy  the  Bragg  matching  condition,  the  vector  triplet  ( Kg,kd,ki )  should 
form  a  closed  triangle  such  that:  kd  =  Kg  +  h.  When  this  occurs  £  will  be  iden¬ 
tically  0.  The  parentheses  in  eq.  3.10  are  intended  to  suggest  grouping  the  differ¬ 
ence  vectors  between  the  object  and  diffracted  wavevector  and  also  the  incident 
and  recording  reference  wavevector.  When  the  difference  vectors  are  identically 
zero,  there  is  strong  diffraction;  hence  the  ideal  reconstruction  of  the  original 
object  beam.  As  a  difference  accrues  between  the  wavevectors,  the  diffracted 
amplitude  decreases.  To  understand  conditions  placed  on  the  mismatch  vector, 
we  continue  the  simplification  of  eq.  3.9  by  making  the  typical  assumption  that 
the  transverse  dimensions  (width  and  height)  of  the  volume  holographic  medium 
are  much  greater  than  the  spatial  bandwidth  of  the  object  beam.  In  this  case 
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Figure  3.1:  Crystal  geometry  and  coordinate  system  for  holographic  recording 
in  the  90  degree  configuration.  An  example  shows  the  superposition  of  scattered 
wavelets  off  of  the  phase  grating. 


the  integrals  over  the  transverse  dimensions  collapse  into  delta  functions  on  the 
— > 

components  of  £  in  the  x  and  y  directions.  The  z  direction  remains  producing: 


A(ki,kd )  oc  f  e  l^'Tdf 
Jv 

poo  poo  pL 

J  —  oo  J  — oc  J  z—0 

=  [‘ 

Jo 


e  1*Tdxdydz 


e~*-*dz 


=  sinc(f*L/2) 


(311) 

(3.12) 

(3.13) 
(314) 


where  sinc(v>)  is  defined  as  sinc(v’)  = 
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3.2.2  Angular  Selectivity 

As  an  example,  let  the  recording  reference  beam  propagate  at  an  angle  9  with 
respect  to  the  x  axis.  The  reconstruction  beam  is  at  a  slightly  different  angle: 
0  +  AO.  The  Bragg  mismatch  in  the  z  direction  is  thus: 

£z  =  ksin(0  +  AO)  -  ksinO  (3.15) 

=  kAOcosO  (3.16) 


where  the  last  step  assumes  A6  -C  9  so  that  a  first  order  Taylor  series  expansion 
of  sin(0  +  A 0)  around  0  is  valid.  The  amplitude  of  the  diffracted  field  in  the 
direction  of  the  object  wavevector  for  a  reference  beam  detuned  by  a  small  angle, 
AO  is: 


A(ki,  kd) 


=  sinc(f2L/2) 

(3.17) 

■  <kML  m 

=  sinc( — - — cos  0) 

(3.18) 

=  s™c(Ae) 

(3.19) 

where  A0  is  the  null  spacing  of  the  sine  function.  Note  that  in  terms  of  the 
wavelength,  A0  is: 


A© 


A 

LcosO 


(3.20) 


If  A  =  0.5  ^m,  L  =  1  cm,  0  =  0  then  A©  =  5  x  10-5  radians  or  2.86  x  10~3 
degrees  (measured  externally  to  the  crystal).  Fig.  3.2  plots  the  angular  selectivity 
for  this  example. 
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Reference  beam  angular  detuning  [degrees] 

Figure  3.2:  Bragg  angular  selectivity  for  a  1  cm  long  crystal  with  A  =  0.5  /xm 
and  a  recording  reference  beam  angle  of  0  degrees. 


3.3  Bragg  Selectivity  Function 


The  general  extension  of  the  angular  selectivity  is  the  Bragg  selectivity  func¬ 
tion  (BSF).  It  yields  the  diffraction  efficiency  as  the  reference  beam  is  detuned 
away  from  the  value  used  during  recording.  The  detuning  may  occur  due  to 
changes  in  wavelength  (A-mux)  or  rotations  of  the  medium  (peristrophic-mux), 
or  a  variety  of  others  possibilities.  Regardless,  in  every  case  it  is  the  change  in 
the  Bragg  momentum  vector,  £,  that  matters. 

In  the  case  of  angular  multiplexing,  we  saw  that  the  transverse  components 
of  the  mismatch  were  constrained  to  be  zero.  In  this  instance,  the  BSF  is  re¬ 
duced  to  a  one-dimensional  function  of  £z,  the  momentum  mismatch  nominally 
in  the  direction  of  the  object  wavevector.  The  diffraction  efficiency  is  thus  the 
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BSF  evaluated  at  f2.  The  strict  relationship  between  £2  and  the  reference  beam 
angle  depends  on  the  specific  interacting  object-reference  grating  vector,  but  the 
dependence  is  slight.  Writing  the  relation  as  in  eq.  3.16,  is  accurate  to  third  order 
for  small  deviations  from  6. 

3.3.1  Material  Effects  and  the  BSF 

The  ideal  case  assumed  the  spatial  modulation  envelope  of  the  permittivity 
to  be  constant.  In  practice,  this  does  not  occur  due  to  the  presence  of  material 
absorption  among  other  things. 

In  this  section  we  wish  to  incorporate  the  attenuated  power  profiles  of  the 
object  and  reference  beams  during  recording  and  also  on  reconstruction.  We  also 
wish  to  assume  an  apodized  form  for  the  reference  plane  wave  as  described  next. 

During  reconstruction  we  choose  R(z)  to  be  the  apodized  field  profile  of  the 
reference  beam.  It  is  defined  such  that  max  |R(z)|  =  1,  thus  the  peak  magnitude 

Z 

of  the  wave  is  Ro.  During  recording,  the  reference  beam  is  defined  as: 

Er(f)  =  W0W{z)eii‘rT~  (3.21) 

where  a  small  reference  beam  angle  is  assumed  in  order  to  ignore  the  geometric 
shadowing  effect  in  W{z).  This  is  valid  to  third  order  deviations  in  the  reference 
beam  angle  away  zero  degrees. 

The  reference  beam  interferes  with  each  plane  wave  component  in  the  object 
beam.  We  model  each  component  in  electric  field  as: 


ik0’f 


E0(f)  =  S0e‘ 


(3.22) 


83 


As  is  common  in  the  literature,  we  assume  all  inferred  angles,  wavelengths,  and 
wavevectors  are  specified  inside  the  material.  We  tacitly  assume  that  the  space 
surrounding  the  hologram  is  index  matched  to  the  crystal.  If  external  angles  are 
desired,  then  Snell’s  law  must  be  applied. 

On  reconstruction,  we  use  a  new  reference  beam  (the  same  as  the  recording 
beam  for  Bragg-matched  conditions)  of: 

Ei(x  =  0,y,z)  =  RoR(z)eikiF  (3.23) 

The  reference  beam  diffracts  through  the  volume  hologram  for  x  €  [0,  W)  and 
2  £  [0,  L}.  As  it  propagates  nominally  in  the  x  direction,  it  will  be  attenuated 
due  to  the  crystal  absorption.  The  reference  beam  at  a  point  t  in  the  medium  is 
expressed  as: 

Ei{f)  =  R0R(z)e-?xeikif  (3.24) 

where  a  is  the  irradiance  bulk  absorption  coefficient.  The  same  relationship  holds 
for  the  recording  reference  beam,  the  object  wave,  and  the  diffracted  wave.  Note 
we  will  assume  only  phase  holograms  so  there  is  no  spatial  modulation  of  the 
absorption  coefficient. 

The  permittivity  modulation  of  the  hologram  results  indirectly  from  the  space 
charge  field  established  during  the  holographic  exposure.  The  fast  fringe  patterns 
in  irradiance  are  mimicked  by  phase-shifted  spatial  electric  fields  established  by 
electrons  bound  in  doped  deep  impurity  traps.  The  linear  electro-optic  effect 
then  produces  the  permittivity  modulation  from  the  space  charge  field.44  For  a 
small  amplitude  modulation  of  the  permittivity,  the  material  index  of  refraction 
approximately  duplicates  the  spatial  distribution.  For  the  90  degree  angular 
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multiplexing  configuration  in  LiNb03,  the  c-axis  is  oriented  at  45  degrees  to  the 
x-z  plane  to  maximize  the  projection  of  the  grating  vector  on  the  electro-optic 
coefficient,  ri3.  In  this  case,  the  permittivity  modulation  is: 

e(r)eiR°f  =  -eo<r13£i  (r)  (3.25) 

where  n0  is  the  ordinary  index  of  refraction  (around  2.286  for  A  =  514  nm)  and 
Ei  ( r)  is  the  space  charge  field  generated  during  the  recording  process. 

Combining  eqs.  3.2  and  3.25  we  compute  the  amplitude  of  the  diffracted  plane 
wave  component  traveling  with  the  correct  object  wavevector,  k0.  The  deviation 
of  the  reconstruction  reference  beam  from  the  recording  beam  is  captured  in  the 
momentum  mismatch  term  £: 

Ed{i)  =  [  (3.26) 

Jv 

=  -e0ntr13  [  RGR(z)e-^L-z+x)ei^-^fE1{rf)dr  (3.27) 

Jv 

=  -eo n\rn  [  RQR{z)e~^L-z+x)ei^-^)fEi{r)e-i^dr  (3.28) 

Jv 

We  see  that  the  diffracted  amplitude  can  be  represented  as  a  three-dimensional 
Fourier  transform  with  f  and  £  forming  a  reciprocal  vector  pair.  Ignoring  some 
of  the  normalization  constants  we  can  write  the  relationship  as: 

Ed(t;)  oc  FT  ^R0R(z)e-^L-z+x)ei^-^)  fE1{r)j  (3.29) 

Remembering  for  the  case  of  angular  multiplexing  that  the  transverse  mis¬ 
match  components  £x  and  must  be  identically  zero,  we  can  simplify  the  rela¬ 
tionship  into  a  one-dimensional  Fourier  transform  pair: 

pH  pW 

Ed(Zz)  oc  /  /  e-^xei{kr~ko)  iEi{r  -x)dxdy  ■ 

Jo  Jo 

FT  ^R0R(z)e-^L-z)ei^r-i‘o)-iEi(r-  z) 


(3.30) 
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We  define  the  Bragg  selectivity  function,  and  the  Bragg  profile  function 
(BPF),  T(z)  from  eq.  3.30  as: 


(3.31) 

(3.32) 


The  BSF  and  BPF  will  be  simplified  once  we  have  a  more  specific  form  for 
the  space  charge  field,  E\{r).  The  first  step  in  finding  E\{r)  is  to  compute  the 
modulation  depth  established  during  recording. 


3.3.2  Modulation  depth 

The  space  charge  field  will  be  proportional  to  the  interference  grating  formed 
between  the  apodized  reference  beam  (defined  in  eq.  3.21)  and  one  of  the  plane 
wave  components  present  in  the  object  beam  (eq.  3.22). 

Interfering  two  plane  waves  with  wavevectors  kQ  and  kr  produces  an  interfer¬ 
ence  pattern  with  a  grating  vector  of: 


(3.33) 
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The  local  irradiance  at  a  location  f  in  the  crystal  is: 


J(r) 

oc  |  E0(f)e~%z  +  Er(r)e-%x\2 

(3.34) 

=  I0(f)  +  E0{r)E*{r)e~^z+3F>  +  Er(f)E*0(r)e~  f(2+I> 

(3.35) 

=  /„(0(l+"(^l('+I,e^+ 

m  ) 

(3.36) 

=  70(f)  (l  +  m(T)etRa'f  +  m*(r)e-li?s'r) 

(3.37) 

m(f) 

S0W0W{z)e-^z^ 

h{r) 

(3.38) 

where  70(r)  is  defined  as: 

Jo(*0  =  \E0^\2e-az  +  \Er(r)\2e-ax  (3.39) 

=  S%e-az  +  WZ\W(z)\2e~ax  (3.40) 

For  short  exposure  times  the  space  charge  field  amplitude  will  grow  linearly 
with  the  modulation  depth  m(r): 

Ei(r;  t)  oc  m{r)eiR°  Ft  (3.41) 

Notice  from  eq.  3.37  that  there  are  possibly  two  interference  gratings  formed. 
One  that  is  Bragg-matched  to  reconstruction  with  the  recording  beam  (ra(r)) 
and  a  second  matched  to  reconstruction  with  the  phase  conjugate  of  the  record¬ 
ing  wave  (m*(f)).  In  a  thick  hologram,  the  conjugate  wave  will  not  diffract  when 
reconstructing  with  a  reference  beam  propagating  in  the  approximate  same  di¬ 
rection  as  the  original  reference  beam,  allowing  us  to  ignore  the  m*(r)  grating 


term. 
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3.3.3  Bragg  Profile  Function 


From  the  modulation  depth  given  in  eqs.  3.36  and  3.37,  and  the  grating  am¬ 
plitude  (eq.  3.41),  we  can  simplify  the  Bragg  profile  function: 

pH  pW 

T(z)  oc  /  /  R0R(z)e-^L-z+x'>ei^-^rE1(r)dxdy 

Jo  Jo 

pH  pW 

=  R0R{z)e-^L-z+^e-i^fm(r)ei^-fdxdy 

Jo  Jo 

pw 

=  H  RoR{z)e~^('L~z+x^m(x,z)dx 
Jo 


S0RoW0R{z)W(z)e-axe-iL 


-dx 


p\V 

H  Jo  S$e~az  +  W2\W(z)\2e~<*x 

H(So\  [w  RoR{z)W{z)e~axe-iL 
\W0)  Jo  (S0/W0)2e-°*  +  \W(z)\*e-a* 

=  H  (jj|)  RoR(z)W(z)e-iL 

pW 

Jo  i 


dx 


(3.42) 

(3.43) 

(3.44) 

(3.45) 

(3.46) 


?  ax  dx 


o  {So/Wo)2e~az  +  \  W{z)\2e~ax 
RoR(z)e~  fr 

I^WI 

IEWli+  (So/»o)V 


(3.47) 


-  H  (  S<1\ 

~  Wo) 


W  +  -  In 


a  ”  \W(z)\2  +  (S0/W0)2e~azeaW 

1  1  +  A(z)e' 

W  +  -  In  w 


(S<L) 

R0R(z)e  *L 

\Wo) 

1  i^wi  L 

a  l  +  A{z)e-a(*-w)  \ 


A(z)  ±  {So/Wo)2/\W{z)\2 

The  integral  in  eq.  3.47  can  be  solved  in  closed  form  as  long  as  W ( z )  ^  0: 


(3.48) 

(3.49) 

(3.50) 


l 


W  e-ax 


-dx  =  [aW  +  ln(C  +  D)  -  ln(C  +  DeaW)]  (3.51) 


I  o  Ce~ax  +  D  aC 

It  is  clear  that  when  W{z)  =  0  the  integral  simplifies  to 

e 


pW 

Jo  I 


~ai  dx  1  -  e~aW 


o  ( So/Wo)2e~*z  +  |IF(z)|2e-“*  (S0/W0yae- 


(3.52) 
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A(z)  is  the  local  apodized  beam  ratio.  The  BPF  and  BSF  defined  in  eqs.  3.31 
and  3.32  can  now  be  expressed  in  their  full  form: 

rL  /  c  N  RoR{z ) 


f  L(So\ 
,(t)  “  Jo  WJ 
=  ft  [r(z>] , 


W  +  1  in  1+A(2)e' 


\W(z)\  L  «  1  +  A(z)e~0‘(z~w'> \ 


e~^zdz  (3.53) 
(3.54) 


where  T(z)  is  defined  to  be  zero  outside  the  interval  [0,  L]  and  also  for  any  2  such 
that  W(z)  =  0. 

We  refer  to  r(^)  as  the  Bragg  profile  function  (BPF)  because  it  is  the  Fourier 
reciprocal  function  to  the  BSF.  Because  the  BPF  is  limited  to  a  spatial  extent  of 
L,  we  see  an  immediate  consequence  is  that  the  BSF  can  not  provide  a  narrower 
main  lobe  than  £0  =  27 r/L.  This  can  be  made  clearer  by  examining  the  Fourier 
transform  pair: 


r(*)  =  1  (3.55) 

vitz)  =  sinc(£zL/2)  (3.56) 

£0  represents  the  smallest  null  spacing  possible  for  the  given  length  L.  It 
thus  establishes  a  practical  minimum  separation  between  successive  holograms. 
Comparing  this  result  to  sec.  3.2.1  we  see  that  this  case  occurs  when  the  overall 
system  acts  like  an  ideal  thick  phase  hologram. 

To  provide  some  real  numbers,  consider  angular  multiplexing  in  a  1  x  1  x  1 
cm3  crystal.  The  ideal  BSF  for  the  9  =  0  page  has  periodic  nulls  at  integer 
multiples  of  Ak  =  sin-1  (A /L).  This  represents  a  nearly  constant  angular  spacing 
of  X/L  «  5  x  10-5  radians  or  2.86  x  10-3  degrees.  There  can  be  up  to  14,000 
multiplexed  holograms  in  a  ±20  degree  span.  Of  course,  such  a  large  angular 
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range  places  additional  non-trivial  design  requirements  on  the  optical  system 
implementation  which  we  tacitly  assume  have  been  adequately  satisfied  (in  this 
chapter  at  least). 

3.4  Apodization  Theory 

The  periodic  nulls  of  the  BSF  are  not  actually  nulls  in  reality  due  to  the 
presence  of  absorption.  Absorption  can  not  be  eliminated  as  it  is  one  of  the  key 
operating  mechanisms  responsible  for  the  photorefractive  effect.  The  first  few 
absorption-corrupted  nulls  can  often  produce  an  unacceptably  large  diffraction 
efficiency  leading  to  significant  interpage  crosstalk.  In  practice  it  is  common  to 
reduce  the  crosstalk  by  storing  pages  not  on  the  first  nulls  of  the  BSF  but  on  the 
second,  third,  or  fourth.  This  solution  markedly  reduces  the  multiplexing  density. 

Another  technique  that  has  been  proposed  to  reduce  the  crosstalk  without 
sacrificing  the  density  is  apodization.13’17’45-47  The  reference  beam  is  shaped 
during  recording  and  reconstruction  in  order  to  provide  a  suitably  improved  BSF. 
Apodization  can  be  applied  in  three  ways: 

1.  during  recording  only:  R(z)  =  1 

2.  during  reconstruction  only:  W(z)  =  1 

3.  during  both  recording  and  reconstruction 

Figure  3.3  shows  the  apodization  profiles  for  (a)  no  apodization,  (b)  recon¬ 
struction  only,  and  (c)  recording  only,  where  the  apodizers  are  chosen  to  flatten 
out  their  associated  r(z).  Figure  3.4  shows  f(z)  for  the  three  apodizers.  For  these 


90 


Figure  3.3:  Apodizer  profiles  for  (a)  no  apodization  (b)  reconstruction  beam 
apodization  (c)  recording  beam  apodization. 


figures  the  assumed  beam  ratio  is  1,  L  =  W  =  H  =  1  cm,  and  a  =  0.73  cm-1. 
Notice  that  the  reconstruction  apodizer  merely  corrects  for  the  nearly  constant 
slope  present  in  the  no  apodization  case.  The  recording  apodizer  accounts  for  the 
exponential  absorption  of  the  signal  beam  as  it  propagates  in  the  +z  direction. 

The  BSF  for  the  three  apodizers  are  shown  in  figure  3.5.  The  non-flat  profile 
of  the  no  apodization  T(z)  results  in  a  relative  amplitude  diffraction  efficiency 
of  5.817%  and  2.894%  at  the  first  and  second  nulls,  respectively.  Reconstruction 
apodization  corrects  most  of  this  by  flattening  out  T(z).  Note  that  recording 
apodization  can  precisely  flatten  T(z)  and  recover  the  perfect  Bragg  nulls. 

During  reconstruction,  all  the  pages  in  the  memory  contribute  a  small  diffracted 
component  superimposed  with  the  Bragg-matched  page.  This  coherent  super¬ 
position  constitutes  a  data-dependent  coherent  field  noise  source  that  we  refer 
to  as  interpage  crosstalk.  When  no  apodization  is  performed  the  nearby  few 
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(a)  (b)  (c) 


Figure  3.4:  Effective  BSF  Fourier  transform  kernel  for  (a)  no  apodization  (b) 
reconstruction  beam  apodization  (c)  recording  beam  apodization. 


Figure  3.5:  Normalized  amplitude  diffraction  efficiency  versus  normalized  multi¬ 
plexing  detuning  for  (a)  no  apodization  (b)  reconstruction  beam  apodization  (c) 
recording  beam  apodization. 
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multiplexed  pages  dominate  the  crosstalk  noise  and  force  a  loss  in  multiplexing 
density  in  order  to  allow  reliable  page  retrieval.  The  above  proposed  recording 
or  reconstruction  apodizers  strongly  reduce  the  cross  talk  noise,  but  there  is  an 
associated  loss  in  diffraction  efficiency  due  to  apodization;  apodization  lowers  the 
system  M/#  as  is  discussed  later  in  sec.  3.6. 


3.4;  1  Experiment 

To  verify  that  we  could  modify  the  Bragg  selectivity  function  by  changing 
the  reference  beam  profile,  we  conducted  some  experiments.  The  optical  setup 
was  constructed  in  a  90  degree  angular-multiplexing  configuration.  A  Newport 
rotation  stage  (0.001  deg  accuracy)  was  used  for  coarse  angular  selection  as  well 
as  for  angular  scans  at  an  approximate  constant  velocity.  A  New  Focus  picomirror 
mount  (0.1  mdeg  accuracy)  was  used  for  fine  angular  control. 

The  photorefractive  crystal  was  placed  3  mm  behind  the  Fourier  plane  of  the 
object  arm,  resulting  in  Fresnel  holograms.  For  this  experiment,  a  5  mm  x  5  mm 
x  5mm  KNSBN  crystal  was  used.  The  crystal  was  chosen  for  its  fast  recording 
time  and  low  absorption  (a  =  0.55  cm-1). 

Figs.  3.6  and  3.7  show  the  reference  beam  profile,  hologram  modulation  depth, 
and  Bragg  selectivity  function  for  rectangular  and  triangular  apodization.  The 
non-uniform  (decaying  exponential)  modulation  depth  in  the  rectangular  apodiza¬ 
tion  clearly  shows  the  effect  of  the  crystal  absorption.  The  “filling-in”  of  the  Bragg 
nulls  as  a  result  is  likewise  apparent  in  fig.  3.6c.  For  the  case  of  triangle  apodiza¬ 
tion,  there  are  no  solid  conclusions.  The  periodic  Bragg  nulls  can  not  clearly  be 


Normalized  power  diffraction  efficiency 


(c)  Bragg  selectivity  function 


Figure  3.6:  Rectangle  apodization.  Experimental  (a)  reference  beam  profile,  (b) 
resulting  modulation  depth,  and  (c)  Bragg  selectivity  function  in  a  5mm  sample 
of  KNSBN.  The  recording  beam  ratio  was  10:1. 


(a)  Reference  beam  profile 


(b)  Modulation  depth 


External  angle  [mdeg] 


(c)  Bragg  selectivity  function 


Figure  3.7:  Triangle  apodization.  Experimental  (a)  reference  beam  profile,  (b) 
resulting  modulation  depth,  and  (c)  Bragg  selectivity  function  in  a  5mm  sample 
of  KNSBN.  The  recording  beam  ratio  was  10:1. 


Irradiance  [a.u.] 
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(a)  Reference  beam  profile 


(b)  Modulation  depth 


(c)  Bragg  selectivity  function 


Figure  3.8:  Gaussian  apodization.  Experimental  (a)  reference  beam  profile,  (b) 
resulting  modulation  depth,  and  (c)  Bragg  selectivity  function  in  a  1cm  sample 
of  Fe:LiNb03.  The  recording  beam  ratio  was  10:1. 
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distinguished.  Though,  we  do  see  the  profile  and  modulation  depth  behave  as  we 
expect. 

Fig.  3.8  shows  the  results  of  using  no  intentional  apodization  (implicitly  a 
Gaussian  apodizer).  The  available  high  beam  power  allows  very  good  resolution 
on  the  sidelobes  of  the  BSF.  On  top  of  the  experimental  curve  is  plotted  the 
theoretical  prediction  for  the  1  cm  Fe:LiNb03  crystal  with  absorption  coefficient 
3.8  cm-1. 


3.4.2  Crosstalk  noise  sources 

So  far  we  have  discussed  apodization  with  respect  to  the  effects  of  the  crystal 
absorption  during  the  recording  and  readout  processes.  Even  when  the  BSF  is 
the  ideal  sine  function,  there  will  still  be  an  intrinsic  amount  of  cross  talk  that 
grows  with  the  number  of  pages  in  the  memory.  It  is  a  directly  a  result  of  the 
non-zero  angular  span  of  object  wavevectors  in  the  x  direction,  as  well  as  a  small 
change  in  the  null  spacing  with  increasing  off-axis  reference  beam  angles.  For 
small  object  spatial  bandwidths  and  reference  angular  range,  this  crosstalk  can 
be  negligibly  small.  But  when  the  number  of  multiplexed  holograms  is  very  large, 
then  the  additive  field  noise  will  offer  a  significant  contribution. 

In  addition  to  intrinsic  cross  talk  noise,  there  are  primarily  two  other  mecha¬ 
nisms  by  which  cross  talk  can  limit  the  system  performance.  First,  if  absorption 
is  not  corrected,  the  neighboring  holograms  will  cross  talk  significantly.  As  dis¬ 
cussed  earlier,  besides  apodization,  this  effect  can  be  reduced  by  spacing  the 
holograms  on  a  larger  angular  spacing,  but  this  may  not  be  a  good  system  level 
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choice  due  to  the  extra  reference  beam  angular  bandwidth  that  is  sacrificed.  The 

larger  angles  also  increase  the  effect  of  the  intrinsic  crosstalk. 

A  second  cross  talk  source  is  a  consequence  of  non-ideal  angular  pointing 

accuracy  in  the  reference  beam.  Errors  in  the  reference  beam  angle  will  decrease 

the  diffraction  efficiency  of  the  desired  page  and  increase  the  cross  talk  from  the 

other  holograms.  The  cross  talk  is  very  sensitive  to  angular  beam  jitter,  and 

therefore  very  small  pointing  errors  will  quickly  magnify  the  cross  talk  levels. 

We  propose  apodization  as  a  technique  that: 

1.  compensates  for  the  non-uniform  grating  profile  (and  hence  restore  the  sine 
nulls) 

2.  increases  tolerance  to  angular  jitter  by  decreasing  the  slope  of  the  BSF 
around  the  nulls 

We  measure  the  effectiveness  of  apodization  by  introducing  the  noise-to-signal 
ratio  (NSR)  which  characterizes  the  cross  talk  noise  by  the  first  and  second  mo¬ 
ments  on  the  detection  plane  irradiance  signal.  We  discuss  a  variety  of  apodiza¬ 
tion  shapes  to  reduce  the  NSR  in  the  presence  of  intrinsic  cross  talk  noise,  ab¬ 
sorption,  jitter,  and  both  a  coherent  and  incoherent  page-to-page  phase  model. 
Afterwards  in  sec.  3.6,  we  derive  the  relationship  between  apodization  and  the 
holographic  system  metric  M /#. 

3.5  Noise-to-Signal  Ratio 

The  diffracted  field  present  at  the  detection  plane  can  be  written  as  a  super¬ 
position  of  every  stored  hologram  weighted  by  their  diffraction  efficiency.  For 
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retrieval  of  page  j  we  write  the  field  amplitude  as: 

Mi 

E  =  Y,  (3.57) 

m=—M  i 

where  the  shorthand  notation  of  £(m;j)  refers  to  the  component  of  the  Bragg 
mismatch  vector  in  the  z  direction  between  pages  m  and  j.  The  large  transverse 
dimensions  (with  respect  to  the  spatial  bandwidth  of  the  object  beam)  require  the 
other  components  of  £  to  be  identically  zero.  This  introduces  a  small  shift  in  the 
diffracted  images  that  we  can  ignore;  however  it  also  imparts  a  small  deviation 
to  £z  as  well,  which  is  accounted  for  in  our  simulation  model. 

The  last  term  in  eq.  3.57  may  or  may  not  be  present  depending  on  the  model. 
The  term  e*^m  represents  a  random  phase  present  upon  retrieval  of  a  page.  When 
we  incorporate  this  term,  we  will  refer  to  the  cross  talk  as  incoherent  (as  the 
average  effect  of  the  term  is  to  average  out  the  cross  talk) .  The  phase  is  assumed 
to  be  uniform  over  the  full  range  of  2 ir.  When  the  term  is  not  present,  the  cross 
talk  is  said  to  be  coherent. 

For  Fourier  holograms,  we  can  write  the  object  beam  as  a  collection  of  plane 
waves  with  direction  cosines  corresponding  to  the  associated  pixel  location: 

E(o]  =  amlpe^pf  (3.58) 

kip  =  A; (cos  ap  cos  fy  z  -  sin  ap  cos  /3ty-  sin  #  x)  (3.59) 

where  ap  and  (3i  are  the  direction  angles  between  the  pixel  at  location  (l,p)  on 
the  SLM  and  the  center  of  the  Fourier  transform  lens  a  distance  /  away  from  the 
SLM.  amlp  is  the  SLM  pixel  value  for  page  m,  pixel  location  (l,p).  In  this  chapter 
we  will  only  consider  a  binary  amplitude  SLM  so  amip  is  either  0  or  1. 
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The  Bragg  mismatch  vector  associated  with  reconstructing  the  grating  cor¬ 
responding  to  page  m  and  pixel  (l,p)  (kr  =  Bm\k0  =  klp),  using  the  reference 
beam  for  page  j  (h  —  k^)  is: 

f  =  k0-kd  +  ki-kr  (3.60) 

=  kip  -kd  +  £(j)  -  k (m)  (3.61) 


where  the  requirements  for  non-zero  diffraction  imply  that  £x  =  0  and  £y  =  0. 
In  order  to  achieve  this,  a  slight  rotation  of  kd  is  necessary  which  imparts  a 
small  change  to  £z.  This  slight  dependence  on  the  x  component  of  the  object 
wavevector,  klp  is  responsible  for  different  columns  on  the  same  page  having 
slightly  different  mismatches  and  hence  slightly  different  cross  talk  values.  We 
can  ignore  the  y  dependence  of  the  mismatch  as  it  is  in  the  degenerate  direction 
(with  an  exception  for  fractal  multiplexing).  As  a  result  of  this  degeneracy,  all 
the  pixels  in  the  same  column  on  the  same  object  page  will  have  similar  cross 
talk  properties.  The  amplitude  of  each  individual  cross  talk  contribution  is  just 
the  Bragg  selectivity  function  evaluated  at  £z. 

The  important  thing  to  note  is  that  the  cross  talk  properties  vary  across 
columns  of  pixels  in  the  detection  plane  as  well  as  from  page  to  page.  Previous 
studies13,46,48  have  analyzed  the  cross  talk  over  all  pages  and  all  columns  and 
found  that  the  worst  cross  talk  occurs  at  approximately  m  =  Mr  =  ±0.9Mi  for 
the  column  of  pixels  with  the  largest  wavevector  component  in  the  x  direction. 
We  will  refer  to  all  cross  talk  measurements  with  respect  to  this  reference  location 
of  page  Mr.  The  diffracted  field,  E,  at  a  pixel  in  the  detection  plane  for  the  worst 
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case  cross  talk  position  can  be  written  as: 

Mi 

E  «  £  (WjKJe*"  (3.62) 

m—-M\ 

6  =  +  (3-63) 

where  (/,  p)  corresponds  to  the  corner  pixel  on  the  SLM. 

Because  of  the  requirement  that  £x  =  0  and  £y  =  0,  there  will  only  be  one 
pixel  grating  from  each  page  that  contributes  to  the  reconstructed  field  for  the 
page  and  pixel  under  consideration.  We  can  simplify  the  notation,  by  implicitly 
assuming  all  the  other  necessary  indices  and  vector  computations,  only  leaving 
the  page  index: 

E  oc  Y2  amT}mel4,m  (3.64) 

m 

By  separating  out  page  MT  from  the  sum,  we  can  write  the  output  field  as  the 
desired  data  embedded  in  an  additive  field  noise: 

E  oc  aMr  +  (3-65) 

m^Mr 

=  aMr  +  En  (3.66) 

Because  the  data  bits  are  random,  the  field  noise  is  likewise  random.  The  first 

and  second  moments  will  strongly  define  how  easy  or  difficult  it  is  to  distinguish 
the  case  of  aMr  =  1  from  aMr  =  0.  We  characterize  this  by  defining  the  noise-to- 
signal  ratio  as: 

NSR  =  (3.67) 

Mi  -  Mo 

where  Hj  and  <jJ  are  the  first  and  second  moments  of  the  detected  pixel  irradiance. 
j  refers  to  the  stored  0  or  1  pixel  value. 


Figure  3.9:  NSR  in  upper  right  quadrant  of  m  =  MT  page. 


The  pixel  irradiance  is  proportional  to  the  square  of  the  electric  field: 

I  oc  \E\2  (3.68) 

=  ^Mr  +  l-En|2  +  2oMrRe  [i?n]  (3.69) 

Fig.  3.9  shows  a  simulation  of  the  NSR  with  respect  to  every  pixel  in  the  upper 
right  quadrant  of  the  m  =  Mr  page.  Notice  that  columns  have  the  same  NSR 
value  as  we  expect  due  to  the  degeneracy  in  the  y  direction. 

3.5.1  Apodizer  Shapes 

We  considered  a  variety  of  apodized  reference  beam  shapes  in  order  to  char¬ 
acterize  the  effect  of  many  parameters  relevant  to  apodization,  such  as  sidelobe 
suppression,  main  lobe  width,  and  sidelobe  decay  rate. 
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In  this  study  we  only  consider  amplitude  apodization  functions.  Phase  apodiza- 
tion  is  a  possible  alternative  technique  which  we  did  not  investigate.  Table  3.1 
lists  the  considered  apodizer  shapes  which  are  described  by  their  normalized  am¬ 
plitude  transmittance  function,  t(z),  defined  over  the  interval  [-5, +5].  Note 
that  some  of  the  apodizers  are  parameterized  by  P  E  [0, 1]. 

Figs.  3.10-3.11  plot  the  transmittance  function  and  spectrum  for  the  apodizers 
considered.  The  spectrum  is  shown  as  the  log-magnitude  in  dB.  Notice  that  the 
peak  value  of  the  spectrum  represents  the  lost  power  when  using  an  absorptive 
apodizers  to  construct  the  apodized  electric  field  profile. 


Apodizer  Type 

t(z) 

Rectangle 

1 

M  >  |(1  -  0) 

else 

Trapezoid 

|  1  -  2|z|//9 

Triangle 

l-2\z\ 

Hamming 

0.5435  +  0.4565  cos  2i rz 

Gaussian 

e-\{20z)2 

Raised  Cosine 

cos^  7 rz 

Tukey 

{  1+ cos  ['»'-«] 

\z\  >  P/2 

_li _ _ _ 

else 

Table  3.1:  Apodizer  functions 


3.5.2  NSR  results 

The  noise  statistics  on  the  irradiance  random  variable  will  determine  when 
cross  talk  noise  limits  the  storage  of  additional  pages.  Consider  the  NSR  when  the 
Bragg  selectivity  function  is  determined  by  the  absorption  corrupted  modulation 
depth.  In  fig.  3.12  we  plot  the  NSR  as  a  function  of  the  number  of  pages  in 
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Figure  3.12:  NSR  (random  page  phase)  for  1  cm  crystal  with  a  -  0.55  cm-1. 


memory,  AT  for  a  variety  of  apodizers  in  addition  to  first  null,  second  null,  and 
third  null  sine  storage.  Note  that  all  the  non-rectangular  apodizers  are  implicitly 
absorption-corrected.  We  see  from  the  default  first,  second,  and  third  null  storage 
that  the  NSR  can  be  reduced  by  decreasing  the  multiplexing  density.  But  clearly, 
if  the  NSR  must  be  lowered,  it  is  far  more  powerful  to  correct  the  absorption  and 
use  first  null  storage,  or  use  one  of  the  other  windowing  functions,  such  as  the 
50%  trapezoidal  apodizer. 

When  the  pages  add  coherently  (fig.  3.13),  we  see  the  NSR  grows  much  faster 
than  with  the  random  phase.  This  result  is  expected  as  the  random  phase  term 
serves  to  average  out  much  of  the  cross  talk  contributions. 
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Figure  3.13:  NSR  (coherent  page  phase)  for  1  cm  crystal  with  a  —  0.55  cm 


When  there  is  an  angular  pointing  error  associated  with  achieving  the  record¬ 
ing  and  reconstruction  angles,  then  the  cross  talk  from  neighboring  holograms 
establishes  a  large  noise  floor.  Fig.  3.14  shows  the  histogram  from  an  experiment 
using  a  New  Focus  picomotor  to  repeatably  point  at  the  same  angle  after  random 
movements  away  from  the  desired  angle.  The  standard  deviation  of  the  pointing 
accuracy  was  measured  to  be  7.6  /urad.  This  corresponds  to  3.3%  of  a  Bragg  null 
for  a  0.5  cm  crystal.  The  irradiance  diffraction  efficiency  is  plotted  and  labeled 
on  the  right  axis.  The  histogram  count  is  labeled  on  the  left  axis. 

The  jitter  dominated  IPI  is  shown  in  fig.  3.15  where  the  pointing  error  is  zero 
mean  and  bounded  below  ±3.3%  of  the  null  spacing.  We  see  in  this  case,  the 
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Figure  3.14:  Experimental  jitter  histogram  from  repeated  measurements  on  point¬ 
ing  accuracy  and  repeatability  of  New  Focus  picomotor.  The  standard  deviation 
is  7.6  //rad,  which  corresponds  to  3.3%  of  a  Bragg  null  spacing  for  a  0.5  cm  crys¬ 
tal.  The  quadratic  curve  and  the  right  axis  refer  to  the  Bragg  selectivity  function 
centered  at  the  first  null. 
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Figure  3.15:  Jitter  dominated  IPI  noise;  L  =  0.45  cm,  a  =  0.55  cm  \  3.3% 
angular  jitter. 

cross  talk  does  not  scale  significantly  with  M  as  the  noise  is  saturated  by  the 
relative  diffraction  efficiency  of  the  neighboring  holograms. 

3.6  Apodized  M/# 

In  this  section  we  consider  the  impact  of  absorption  on  the  recording  and  era¬ 
sure  dynamics  of  photorefractive  volume  holograms,  a  is  a  key  component  in 
the  photorefractive  effect  driving  the  time  rate  of  grating  formation  and  erasure. 
When  the  absorption  is  very  small,  holographic  gratings  record  slowly  and  erase 
slowly.  When  the  absorption  is  larger,  gratings  record  faster  but  subsequently 
erase  faster.  Somewhere  between  the  two  extreme  cases  exists  an  optimal  absorp¬ 
tion  coefficient  for  a  given  set  of  material  and  system  parameters.  This  relation 
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was  derived  from  the  underlying  crystal  physics  and  material  parameters  in  ref¬ 
erence  [42],  The  absorption  is  optimal  in  the  sense  that  it  maximizes  the  system 
metric  M/#.31  M/#  provides  the  scaling  constant  relating  the  Bragg-matched 
page  diffraction  efficiency  and  the  number  of  multiplexed  pages  M: 


(3.70) 


where  r)(£z  =  0)  is  the  Bragg-matched  amplitude  diffraction  efficiency  considered 
in  section  3.3.  Following  a  similar  approach  to  reference  [42]  we  develop  the 
M/#  for  the  apodized  system  using  Fe:LiNb03  as  the  photorefractive  medium. 

The  reconstruction  reference  beam  (eq.  3.23)  incident  on  the  crystal  has  a 
total  power,  Pin,  of: 


pH  pi 

Pin  =  /  \RoR(z)\2dydz 

Jo  Jo 

=  HR l  [  \R(z)\2dz 
Jo 


=  HB%Pr 


(3.71) 

(3.72) 

(3.73) 


where  for  convenience  we  define  the  power  in  the  apodizer  shape  as: 


=  /*!*(*)  Is 

Jo 


(3.74) 


The  diffracted  power,  P^t,  can  be  found  by  integrating  over  the  angular  spec¬ 
trum  of  the  diffracted  wavevectors: 

P out  =  //I  A(kd,  ki)sinc(k,itxW)smc(kdtyH) |  dkdtXdkd>y  (3.75) 


no 


which  by  Parseval’s  theorem27  can  be  reexpressed  as  the  power  in  the  reciprocal 
domain: 


Pi out 


nw 

|g«(Q)l2 

WH 


^rect(^)rect  (•§) 


dx  dy 


(3.76) 

(3.77) 


where  Ed{^z)  was  defined  in  eq.  3.28. 

The  amplitude  diffraction  efficiency  is  the  ratio  of  the  square  root  of  the  in¬ 
put  and  diffracted  powers,  allows  us  to  represent  the  Bragg-matched  diffraction 
efficiency  (BSF  at  £2  =  0)  as: 


77(0) 


\EM\ 

HRoVWPR 


€o^r13 


ill 


Ei(r)R(z)Sa(x,z)e 


-iKq 


(3.78) 

(3.79) 

(3.80) 


For  small  recording  exposures  the  space  charge  field  and  the  modulation  depth 
are  related  as: 


| JE7i (r;  i)|  oc  m(x,  z)Esc(a)  (l  - 


e~t/T,(x,z) 


(3.81) 


where  the  new  parameters  are  defined  in  Burr’s  thesis  [42].  r(  is  the  local  recording 
time  constant,  uji  is  the  imaginary  component  of  the  local  recording  time  constant. 
Esc  summarizes  a  number  of  material  properties.  It  can  be  thought  of  as  a 
material  scaling  parameter.  It  is  important  to  note  that  it  depends  on  a  among 
other  things. 
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In  terms  of  the  saturation  amplitude  and  recording  time  constant,  the  space 
charge  field  evolves  according  to: 


|£i(r;*)l  =  Ao[  l-e-*/Tr] 


(3.82) 


Taking  the  time  derivative  at  t  —  0  for  both  forms  allows  us  to  express  the 
recording  slope  ( A0/rr )  as: 


eorc£ri3 


|  JJI  S0W0W{z)R(z)e-^x+Vdr 


(3.83) 


eongr13g0lToe-fL(l  -  e ~aW) 
TxVWPr 


|-£'sc(q:)  (1  +  iTiUJi) | 


f  W(z)R(z)dz 
Jo 


(3.84) 


During  erasure,  the  space  charge  field  decays  as  e-t/Te.  re  can  be  found  by 
taking  the  ratio  of  r]  and  the  derivative  at  t  =  0: 

T  —  Tx  fo  R(Z)W (z)  dz 

5o  To  R(z)W(z)e-°*dz  +  R(z)W{z)\W(z)\*dz 

Finally,  the  M/#  can  be  computed  as  the  product  of  the  recording  slope  and 
the  erasure  time  constant: 


M/#  &  (. A0/rr)re  (3.86) 

=  x  (3.87) 


R(z)W(z)dz j 

2 

1 

( So \ 
ylTo ) 

|  [L  R(z)W(z)e 
Jo 

1  +  e~aVV 

~azdz + 

2 

fL  R(z)W(z)\W{z)\2dz 
Jo 
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where  k0  =  27r/A0,  n0  is  the  ordinary  index  of  refraction,  r13  is  electro-optic 
coefficient,  and  Esc(a )  is  the  space-charge  field  amplitude  which  depends  on  a 
number  of  crystal  parameters  including  a. 

In  figure  3.16  we  plot  the  M/#for  the  three  apodization  cases  versus  the 
absorption  coefficient.  For  each  choice  of  a  the  beam  ratio  is  selected  to  maximize 
the  M/#.  Recording  apodization  restores  the  ideal  Bragg  nulls  but  suffers  a  6.8% 
loss  in  M /#.  Apodization  of  the  reconstruction  beam  loses  1.5%  in  M/#  and 
nearly  perfectly  restores  the  ideal  nulls. 

An  important  result  to  notice  is  that  the  optimum  absorption  coefficient  (that 
maximizes  M/ jfc)  has  shifted  due  to  apodization.  For  a  1  x  1  x  1  cm  LiNb03  crys¬ 
tal,  the  optimal  a  is  0.727,  0.697,  0.608  cm"1  for  no  apodization,  reconstruction 
beam  apodization,  and  recording  beam  apodization,  respectively. 


3.7  Summary 

In  this  chapter  we  have  presented  some  theoretical  results  about  page-to-page 
cross  talk  in  photorefractive  media.  We  saw  the  the  Bragg  selectivity  function 
can  be  related  to  the  transverse  electric  profile  through  a  Fourier  transform.  This 
allows  modification  of  the  reference  beam  to  reduce  cross  talk  and  correct  for 
deleterious  material  effects  such  as  absorption.  In  addition,  we  developed  the 
modifications  necessary  to  describe  the  M/#  of  an  apodized  system.  Simula¬ 
tions  showed  that  there  was  only  a  minor  cost  in  reduction  of  M/ #  when  using 
apodization  to  correct  the  influence  of  absorption  and  restore  the  ideal  Bragg 


nulls. 
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Figure  3.16:  M/#  versus  absorption  in  a  1  x  1  x  1  cm3  Fe:LiNb03  crystal  for  no 
apodization,  reconstruction  beam  apodization,  and  recording  beam  apodization. 

Notes  Portions  of  this  work  have  been  published  in  Proceedings-of-the-SPIE 
as  described  in  reference  [47].  The  basic  derivation  of  M/#  is  due  to  Geoff  Burr.42 
Extensions  to  incorporate  apodization  in  the  derivation  were  originally  considered 
by  a  former  OCPL  group  member,  Xue-Wen  Chen.  There  are  only  a  few  minor 
differences  between  his  and  my  derivation. 


